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We present a formalism for precise description of oscillation phenomena in matter at high energies 
or high densities, V > Am 2 /2E, where V is the matter-induced potential of neutrinos. The accuracy 
of the approximation is determined by the quantity sin 2 26mAV/2irV , where 6 m is the mixing angle 
in matter and AV is a typical change of the potential over the oscillation length (I ~ 2ir/V). We 
derive simple and physically transparent formulas for the oscillation probabilities, which are valid for 
arbitrary matter density profiles. They can be applied to oscillations of high energy (E > 10 GeV) 
accelerator, atmospheric and cosmic neutrinos in the matter of the Earth, substantially simplifying 
numerical calculations and providing an insight into the physics of neutrino oscillations in matter. 
The effect of parametric enhancement of the oscillations of high energy neutrinos is considered. 
Future high statistics experiments can provide an unambiguous evidence for this effect. 



PACS numbers: 14.60.Pq, 14.60.Lm 

1. Introduction. Neutrino physics enters a new phase 
now, where the objectives are precision measurements of 
the parameters, studies of subleading oscillation effects 
and searches for new physics beyond the already stan- 
dard picture, which includes non-zero neutrino masses 
and mixing. Detection of neutrinos from new sources, in 
particular, of cosmic neutrinos, is in the agenda. 

Substantial new information is expected from the stud- 
ies of high energy [E > 1 GeV) neutrinos. This includes 
investigations of atmospheric neutrinos with new large 
volume detectors [[J, long baseline accelerator experi- 
ments and detection of cosmic neutrinos from galactic 
and extragalactic sources, as well as of neutrinos pro- 
duced in interactions of cosmic rays with the solar at- 
mosphere |3- Another possible source of neutrinos is 
annihilation of hypothetical WIMP's in the center of the 
Earth and the Sun |j] . In all these cases beams of high 
energy neutrinos can propagate significant distances in 
the matter of the Earth (or of the Sun) and therefore 
undergo oscillations/conversions in matter. 

Increased accuracy and reach of neutrino experiments 
put forward new and more challenging demands to the 
theoretical description of neutrino oscillations. In the 
present Letter we derive very simple and accurate ana- 
lytic formulas describing oscillations of neutrinos in mat- 
ter. Our primary goal is to study oscillations of high 
energy neutrinos [f| , but the formulas we obtain are ac- 
tually applicable in a wide range of neutrino energies. 
They simplify substantially numerical calculations and 
allow a deep insight into the physics of neutrino conver- 
sions in matter. In particular, they provide a useful tool 



for studying parametric resonance enhancement of neu- 
trino oscillations. 

The parametric resonance can occur in oscillating sys- 
tems with varying parameters due to a specific corre- 
lation between the rate of the change of these parame- 
ters and the values of the parameters themselves. In the 
case of neutrino oscillations in matter, the parametric 
enhancement is realized when the variation of the mat- 
ter density along the neutrino trajectory is in a certain 
way correlated with the change of the oscillation phase 
HQ. This effect is different from the MSW effect 
and yet can result in a strong amplification of the oscil- 
lations. We apply the formalism developed here to study 
the parametric resonance in oscillations of high energy 
neutrinos in the Earth. 

2. Formalism. We consider oscillations in the 3-flavor 
neutrino system (y e , v^, v T ), with the mass squared dif- 
ferences Am^ and Ato^ responsible for the oscilla- 
tions of atmospheric and solar neutrinos, respectively. 
In certain situations, the full three-flavor problem is ap- 
proximately reduced to the effective two-flavor ones, in 
which the electron neutrino v e oscillates into a combi- 
nation of and v r . We shall be mainly interested in 
oscillations of neutrinos with energies E > Am^ 1 /2V, 
where the matter-induced potential of neutrinos V(x) = 
\Z2GfN e (x), with N e (x) the electron number density 
in matter and Gf the Fermi constant. Numerically, 
this corresponds to E > 8 — 10 GeV for the matter of 
the Earth. In this case the 1-2 mixing is strongly sup- 
pressed by matter, and the problem is reduced to an ef- 
fective two-flavor one, described by the potential V , mass 
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squared difference Am 2 = Am^ and vacuum mixing an- 
gle 9 = 6 1 3 (which is assumed to be non-zero) [l(|. In 
particular, the oscillations of electron neutrinos are de- 
termined by the transition probability P 2 = P(y e <-> i/ a ), 
where v a = sin #23^ + cos#23^V is the state which mixes 
with v e in the third mass eigenstate. In terms of P2 the 
flavor transition probabilities are P{y e — > f^) = Piy^ — > 

^e) = Sin 2 023-P2, P(^e ~ * Vt) = P{v T — ► f e ) = COS 2 6>23-P2 

[lCj . The third state (orthogonal to v a ) decouples from 
the rest of the system and evolves independently. 

In the (y e , u a ) basis, the evolution matrix S{x) describ- 
ing neutrino oscillations satisfies the equation 



with the Hamiltonian 



H(x) 



V 
~2 



-cos26» sin 29 
sin 29 cos 29 



Here 



5 = Am 2 /AE , 



(1) 



(2) 



(3) 



and the first (potential) term dominates in the high en- 
ergy limit. However, in most situations of interest the 
neutrino path length in matter L satisfies 5 ■ L > 1, 
therefore we cannot consider the whole second term as 
a small perturbation, and the effect of 5 on the neutrino 
energy level splitting should be taken into account. For 
this reason we split the Hamiltonian as 



with 



Ho 



1 

-1 



H — H + Hi, 



Hi = sin 29 S 



(4) 



(5) 



Here 



u(x) ee sj (V/2 - 6 cos 29) 2 + 5 2 sin 2 29 , (6) 

2u> being the difference of the eigenvalues of H(x); 

cos295-V/2 + uj 8 . nA 1 

6 ee «— sin20<l. (7) 

sin 29 5 V y ' 

The ratio of the second and the first terms in the 
Hamiltonian in (0) is determined by the mixing angle 
in matter m : 



sin 2(9 <5 



sin 29 m ■ 



(8) 



Therefore for sin 28 m -C 1 the term Hi can be considered 
as a perturbation. Furthermore, according to Q, £ ~ 
sin 29 m , so that the diagonal terms in Hi can be neglected 
in the lowest approximation. 

We seek the solution of eq. (QJ in the form 



So ■ Si 



(9) 



where So is the solution of the evolution equation with 
H replaced by H - From |J5J| we find 



Sq(x) 



-i<f>{x) Q 





where 



4>{x) 







dx' u) (x) 



(10) 



(11) 



is the adiabatic phase. Then, according to P|. |1| and 
@, the matrix Si satisfies the equation 



i — ; — — Sn ^HiSo Si — H1S1 
ax 



(12) 



where Hi = S Q 1 HiSq is the perturbation Hamiltonian in 
the "interaction" representation. Eq. (|12|l can be solved 

by iterations: Si = 1 + Sj 1 " 1 + which leads to the stan- 
dard perturbation series for the S matrix. For neutrino 
propagation between x — and x — L we have 

/>L pL />x 

Si{L) = 1-i / dxili{x)- \ dxHi(x) \ dx'Hi(x')... 
Jo Jo Jo 

(13) 

Taking Si to the first order, we obtain from ijjjj the evo- 
lution matrix 



S(L) = S (L) 



1-iS sin 29 / dx 







-i2<p(x) q 



(14) 

The v e «-> v a transition probability Pi (L) is given by the 
squared modulus of the off-diagonal element [S(L)] ae : 



P 2 = S 2 sin 2 29 



dx e ^ 2 ^ x) 



(15) 



For density profiles that are symmetric with respect to 
the center of the neutrino trajectory, V(x) = V(L — x), 
eq. lfT5|l gives 



Pi 



A? 



4E 



sin 2 29 



L/2 



dz cos 2(f>(z) 



, (16) 



where z = x — L/2 is the distance from the midpoint 
of the trajectory and <f)(z) is the phase acquired between 
this midpoint and the point z. 

The transition probability P2 scales with neutrino en- 
ergy essentially as E~ 2 . The accuracy of eq. (fT5|l also 
improves with energy as E~ 2 . This is illustrated by pan- 
els (a) and (b) of fig. 1, which show P2 as a function of 
neutrino energy for several trajectories inside the Earth 
and as a function of the zenith angle of the neutrino tra- 
jectory 6 for several neutrino energies. One can see that 
already for E > 8 MeV the accuracy of our analytic for- 
mula is extremely good. Note that when neutrinos do 
not cross the Earth's core (cos 9 > —0.837) and so ex- 
perience a slowly changing potential V(x), the accuracy 



3 



of the approximation (|15|) is very good even in the MSW 
resonance region E ~5 - 8 GeV. The accuracy of eq. i|15|) 
is also good for energies below ~ 2 GeV (not shown in 
the figure); however, in this region the domain of the ap- 
plicability of (|15|l is relatively narrow, since for E < 0.5 
GeV the oscillations driven by the "solar" parameters 
(Am|i, 612) can no longer be neglected. 

To understand the remarkable accuracy of eq. Q15[l. 
we find the correction AP 2 to the transition probabil- 
ity in (|15|l emerging in the next nontrivial order in Hi. 
Note that from the above considerations one can expect 
AP2/P2 to be proportional to sin26* m . Furthermore, as 
can be easily seen, for uniform matter eq. (|15|l repro- 
duces the exact transition probability; therefore one ex- 
pects AP2/P2 oc V . A straightforward calculation in- 
deed gives 



AP 2 

P2 



sin 2 26V, 



AV 

Attlo 



sin 2 29. „ 



AV 
2nV 



(17) 



where AV is the change of the potential over the oscilla- 
tion length 7r/(jj, and the last equality holds in the high 
energy regime. For slowly changing density this is equiv- 
alent to 



AP 2 



sin 2 20 „ 



V 
Alo 2 



(18) 



Introducing the adiabaticity parameter 7 = 
Attuj /(sm29 m AV), we find that AP 2 /P 2 ^ sm29 m j~\ 
and therefore for small mixing in matter our approxima- 
tion is better than the adiabatic one. At the same time, 
for AV/4irui < 1 it is better than the simple expansion 
in powers of sin 2 29 m . 

The matter density profile of the Earth satisfies 
V /V 2 < 0.5, and therefore for neutrino oscillations in 
the Earth our approximation is expected to work well 
when sin 2 29 m <C 1 . This is fulfilled in the high energy 
(or, equivalently, high density) limit EV/ Am? ^> cos 29, 
i.e. above the MSW resonance. If the vacuum mixing 
angle is small (i.e. 9 = #13), our expansion parameter 
is also small below the resonance. The above formalism 
applies in this low energy case as well, with only minor 
modifications: the sign of Ho in (J5J has to be flipped, 
and correspondingly one has to replace u) — > —lu in eq. 
0. These modifications are necessary because of the 
interchange of the eigenvalues of the Hamiltonian upon 
crossing the MSW resonance. Expressions for the tran- 
sition probability in eqs. (jTBTl . (|H)|l remain unchanged. 
Thus, our results are in general valid outside the MSW 
resonance region, which for small 9 is very narrow. For 
the non-resonant channels [p channels for Am 2 > or v 
channels for Am 2 < 0) and small vacuum mixing our for- 
mulas are valid in the whole diapason of energies because 
sin 29 m is always small. 

If 0i3 is very small or vanishes, v e «-> v^ r oscillations 
are driven by Am 2 = Aniji and the large mixing angle 
9 = 9i2- The oscillation probabilities can then be ex- 
pressed through another effective two-flavor probability, 



P2 = P2(Am 2 1 , 9i 2 , V(x)), in terms of which the fla- 
vor transition probabilities are P(v e <-> v^) — cos 2 ^23^2, 
P(v e <h> v T ) = sin 2 6> 23 P 2 EH- For EVjAm\ x > cos26>i 2 
(which for the typical densities inside the Earth corre- 
sponds to E > 0.5 GeV), the probability P 2 is very well 
approximated by eq. I|15l) . 

Let us consider the case of symmetric matter. Inte- 
grating (|l(jfl by parts, one finds 



P 2 = sin 2 26»° 



sm0i + lu 



L/2 



I 2 



dz — sm2(j)(z) 



duo 1 

dz U! 2 

(19) 

where 9® m = 9 m (V~o), uo = uo(Vq), Vo being the poten- 
tial at the initial and final points of the neutrino tra- 
jectory, and 4>l is the adiabatic phase acquired along 
the entire neutrino path. If the potential changes slowly 
with distance, so that uo~~ 2 duo/dz <C 1, the second term 
in i|T§|) can be neglected, and P 2 reduces to the usual 
adiabatic probability in symmetric matter: P a diab = 
sin 2 29 Q m sin 2 tfii . The second term in (|19|) describes the 
effects of violation of adiabaticity. 

Let us apply the above results to neutrino beams (at- 
mospheric, accelerator, cosmic neutrinos) crossing the 
Earth. According to the PREM model [3, the Earth 
density profile can be described as several spherical shells 
of radii Ri with rather smooth density change within the 
shells and sharp change at the borders between them. 
Then, along a direction from the center of the Earth out- 
wards, uj(z) decreases abruptly from uof to w,r in very 
narrow regions around Ri. Therefore dio/dz is large in 
these narrow regions and small outside them. The inte- 
gration in l|19|l can then be easily done, leading to 



P> = sin 2 26> 







sm <Pl 



ujy — U). 
cjo^i 1 _ l sin 20 4 



(20) 



Here 4>i is the adiabatic phase acquired by neutrinos be- 
tween the points z = and z — Ri. We will use eq. I|20|) 
for interpreting the results of our calculations. 

3. Parametric enhancement of oscillations. In the 
PREM model of the Earth density all the density jumps 
between different shells except those between the mantle 
and core are relatively small E3- Therefore the density 
profile felt by neutrinos crossing the core of the Earth can 
be approximated by three layers (mantle - core - mantle) 
with slow density changes within the layers. Eq. I|20() 
then gives 



P 2 = sin 2 29l 



\{<j) c + 2(j), 



s _ ^0 
1) 



1 



sine 



~] 2 



(21) 



where uj m and u> c are the values of uj(x) in the mantle 
and core on the respective sides of their border, and cf> m 
and 4> c are the phases acquired in the mantle (one layer) 
and core. Eq. I|21() corresponds to the adiabatic neutrino 
propagation inside the core and mantle and violation of 
the adiabaticity at the borders between them. In the ap- 
proximation of constant effective densities in the mantle 
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FIG. 1: Panels (a) and (b): transition probability P2 vs neutrino energy for different trajectories inside the Earth [panel (a)] 
and vs cos O for different neutrino energies [panel (b)] . Solid curves are the results of exact numerical calculations, dashed curves 
are obtained using formula 1151 . Small window (b') in panel (b) shows the values of the parameter X3 [eq. 1251 1 calculated in 
the three- layer model of the Earth's density. Panel (c): contours of constant phases m (solid curves), (j> c (dashed curves) and 
of constant X3 = (dotted curves). The numbers at the curves are the values of the phases in units of tt. The shaded areas in 
panesls (b) and (c) correspond to the Earth's core. For all panels we take sin 20i3 = 0.15 and Am 2 = 2 • 10~ 3 eV 2 . 



w^f , and the non-adiabatic term 

, .off /, ,eff\ 



and core we have ujq 
is proportional to (1 

For neutrino trajectories that cross the mantle only 
((f> c = 0), eq. l(2*Tl) reduces to the adiabatic probabil- 
ity. The passage of neutrinos through the core can lead 
to an enhancement of the oscillations. As follows from 
(|2ip. the maximum enhancement of the probability can 
be achieved when sin(</> c + 20 m ) and sm<p c are of opposite 
sign and maximal amplitude: 



sin 



+ 2</> m ) = ±l, 



that is, when 



7T 7T 
±-+27m, (f) m = ±- 



2nk. 



(22) 



(23) 



Here n and k are integers, and the signs in front of tt/2 
are correlated. In this case the enhancement factor is 



r>max 

;in 2 28„ 



U) 



*2.5, 
(24) 

where sin 2 29 m in the denominator corresponds to the 
maximum possible transition probability for neutrinos 
crossing only the mantle, and we have taken into account 
that at high energies uj m /uj c w V m /V c - 

The condition in eq. I|23|) and the enhancement de- 
scribed by eq. are the particular cases of the para- 
metric resonance condition and the parametric enhance- 
ment of neutrino oscillations 0, • In [13 it was shown 
that in the case of matter consisting of alternating lay- 
ers of two different constant densities and (in general) 



different widths the parametric resonance condition is 



- (sin 4> m cos 4> c cos 28 m + cos < 



where 9 m ,c & n d 



sin (j) c cos2(9 c ) = , 
(25) 

are the mixing angles and the ac- 
quired oscillations phases in the layers m and c. This 
condition can also be used as an approximate one when 
matter density inside the layers is not constant but varies 
sufficiently slowly. For neutrino oscillations in the Earth 
we identify the layers m and c with the mantle and 
core. Since in the energy region sin 28 S <C V one has 
m ps 8 C « 7r/2, condition (|25|) reduces to 



X, 



sin 



=) =0. 



(26) 



Eq. ilL'-'SI) is a particular realization of this condition. In 
the high energy limit the parametric resonance condi- 
tion was previously considered in the framework of 
active-sterile atmospheric neutrino oscillations in |14| . It 
should be noted that, while the realization l|23|) of condi- 
tion (|25|l corresponds to the maximal possible parametric 
enhancement of oscillations of high energy neutrinos in 
the Earth, a sizable amplification is also possible if the 
equality X% — is realized differently, i.e. when the two 
terms in (12511 do not separately vanish but cancel each 
other [l5j. 

The parametric resonance conditions in eq. 1251) or 
(|23fl require a subtle correlation between the matter den- 
sity of the Earth, the distances that neutrinos travel in 
the Earth's mantle and core (which are not independent), 
and also in general neutrino energy and oscillation pa- 
rameters. Therefore it is far from obvious that these 
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conditions can actually be satisfied. Amazingly, this is 
indeed the case . pj . pj . 

Our present analysis shows that, for v e <-> v^ T os- 
cillations of high energy neutrinos in the Earth, there 
are two regions of the zenith angles of neutrino trajecto- 
ries in which (|25(l is satisfied: cosG = — 1 -, — 0.93 and 
cos 6 = -0.88 ~ -0.84. This is illustrated in panel (b) 
of fig. 1. Since at high energies matter suppresses neu- 
trino mixing, one could expect that for the trajectories 
crossing the core (cos 9 < —0.837), where the densities 
are higher, the probability would be suppressed or at 
least would not change. Instead, we see two prominent 
peaks there, exceeding maximal allowed by the MSW ef- 
fect value of probability sin 2 29 m by up to a factor of 
2. This is the result of the parametric enhancement of 
neutrino oscillations. As can be seen from the figure, for 
core-crossing trajectories the positions of the peaks of Pi 
essentially coincide with zeros of A3 . 

For neutrino energies E ~ 10 - 15 GeV, the oscillation 
phases corresponding to the peak with cos 9 < —0.93 
(the inner peak) are <p rn ~ 7r/4, (f> c ~ 77r/4, while for 
the peak with cos 9 = — 0.88 -. — 0.84 (the outer peak), 
they are (j) m ~ 0.357T, <j) c ~ 0.657T. In both peaks to a 
good accuracy (f> m + (f> c — nn, so that eq. l(2l)|> is satisfied. 
The phases in the outer peak are closer to the realization 
(|23[l of the parametric resonance condition, and therefore 
in this peak the parametric enhancement of oscillations 
is closer to the maximal possible one. From panel (c) 
of fig. 1 one can see that at E ~ 21 GeV the maximum 
enhancement condition (|23|l can be exactly realized in the 
outer peak. For neutrinos of very high energies (E > 100 



GeV), it can be realized nearly exactly in the inner peak, 
whereas the outer peak becomes slimmer and lower, and 
at energies E > 300 GeV virtually disappears. 

4- In conclusion - We have derived very simple and ac- 
curate integral formulas describing neutrino oscillations 
in matter with arbitrary density profiles. They can be 
applied to all possible situations where the effective mix- 
ing in matter is small. In particular, our results can be 
applied to atmospheric, accelerator and cosmic neutrinos 
crossing the Earth as well as to cosmic neutrinos cross- 
ing the Sun. They can be used for probing the effects of 
small jumps in the density profile of the Earth as well as 
evaluating the effects of uncertainties in this profile on 
the interpretation of the results of future accelerator ex- 
periments. They can also be employed for studying the 
effects of proper averaging over the energy spectrum of 
the neutrino beam. We used the obtained formulas to 
study the parametric enhancement of oscillations of high 
energy neutrinos in the Earth and identified two peaks in 
the zenith angle distribution of core-crossing neutrinos, 
which are due to the parametric effects. Observation 
of these peaks in future high statistics experiments with 
high energy neutrinos would provide an unambiguous ev- 
idence for the parametric resonance effects in neutrino 
oscillations in matter. 
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